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Abstract 

In this paper we consider a sample of a linearly elastic heterogeneous com- 
posite in elastodynamic equilibrium and present universal theorems which 
provide lower bounds for the total elastic strain energy plus the kinetic en- 
ergy, and the total complementary elastic energy plus the kinetic energy. 
For a general heterogeneous sample which undergoes harmonic motion at a 
single frequency, we show that, among all consistent boundary data which 
produce the same average strain, the uniform-stress boundary data render 
the total elastic strain energy plus the kinetic energy an absolute minimum. 
We also show that, among all consistent boundary data which produce the 
same average momentum in the sample, the uniform velocity boundary data 
render the total complementary elastic energy plus the kinetic energy an 
absolute minimum. We do not assume statistical homogeneity or material 
isotropy in our treatment, although they are not excluded. These univer- 



sal theorems are the dyna mic equivalent of the un i yersa 



known for the static case (iNemat-Nasser and Horil (119951 )). It is envisaged 



theorems already 



that the bounds on the total energy presented in this paper will be used to 
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formulate computable bounds on the overall dynamic properties of linearly 
elastic heterogeneous composites with arbitrary microstructures. 
Keywords: Energy Bounds, Effective Dynamic Properties, Metamaterials, 
Homogenization 



1. Introduction 



One of the main objectives of micromechanics is to estimate the overall 
properties of a heterogeneous composite in terms of the architectural and ma- 
terial properties of its micro- constituents. For the static case, this amounts 
to defining effective properties which relate the domain averages of the stress 
and strain tensors over a suit a bly la r ge samp l e calle d a Repr e sentat i ve Vol 



Willis 



ume E lement (RVE; See IffiU (11963r i: iHashinl (119651 ): iKroned fll977f ) 
( Il981al jbl)). The analysis is complicated by the following twin problems: 

1. It is often difficult to identify a suitable RVE that characterizes the 
composite. 

2. In general the average values of the stress and strain tensors depend 
upon the boundary conditions to which the RVE is subjected. Therefore, ef- 
fective properties defined to relate these averages depend upon the boundary 
conditions under which the averages are calculated. 

To deal with these difficulties, statistical homogeneity is used which as- 
sumes that the overall properties of the composite can be represented by 
those of an RVE. Furthermore, it implies that the overall response of the 
RVE is almost the same for any boundary condition as long as the aver- 
age stress and strain tensors are kept constant. By employing statistical 
homogeneity, strain and complementary energy functionals can be defined 
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to est i mate the overall proper t ies of 



195 



Hashin and Shtrikmanl (Il962a 



;lie heterogeneou s composit e 



Eshelbv 



bh: iHashinl f llDGSh : IWalpolel f ll966al 



Korringal ( 119731 ) . However, such estimates have been shown to be merely 



plausib le bounds and not rigorous bounds on the effective properties 



Willis 



(Il981al ). It was subsequently shown that exact inequalities for the elastic and 



the complementary energies sto red in a finite body under arbitrary b oundary 



conditions could be established 



Nemat-Nasser and Hori 



(11995 



1999|). These 



bounds on the energies were, in turn, used to provide improvable and com- 
putable bounds on the static effective prope rties of heter o geneo us composites. 



The recent interest in metamaterials 



(120031); 



Liu et al 



Smith et al. 



(120051 ): Milton and Wilhsl (120071 ) has necessitated a sys- 



feoOOl); 



Sheng et al 



tematic homogenization procedure to estimate the effective dynamic prop- 
erties of composites by us i ng field averaging or ensemble a veragin g t echniques 



Smith and Pendrvl (120061): 



Nemat-Nasser et al. 



Amirkhizi and Nemat-Nasserl (l2008al lbl) 



Willis 



(120091 ): 



20111 1: IWiUid (1201 iD . It has been shown that in the 
homogenized representation the average stress is coupled with the average 
velocity a nd that the ave r age rn omentum is coupled with the average strain 



(See also 



Shuvalov et al. 



(l201l[ )). Recently, a micromechanical method to 



calculate the effective propert i es of p e riodic elastic composites was p r opose d 



Nemat-Nasser and Srivastaval (120111 ): ISrivastava and Nemat-Nasserl ( 1201 ll ). 



This method provides effective parameters which reduce to the e ffective pa- 



ramet ers calculated from the ensemble averagin g technique of Ref . 



Milton and Willis 



Nemat-Nasser and Srivastava 



(2007) or the field averaging technique of Ref. 
( 120111 ) when the dispersion relation of the composite is enforced. 

In this paper we begin with a brief overview of the elastostatic problem. 
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stating the universal energy theorems for this case. For the elastodynamic 
case, we present the general form of the effective constitutive relations and 
briefly discuss the properties of the effective parameters which arise in the 
averaged constitutive relations. Finally, we present universal theorems which 
are the dynamic analogues of t he static universal theorems presented in Ref. 
Nemat-Nasser and Horil (119991 ). These universal theorems bound the total 
elastic energy plus the kinetic energy and the total complementary energy 
plus the kinetic energy of a composite in an elastodynamic state and pave the 
way for subsequently establishing rigorous bounds on the effective dynamic 
properties of the composite. 

2. Elastostatic Universal Theorems 

Consider the static case of a general heterogeneous solid which consists 
of various elastic phases. There is no restriction on the number, geometry, 
material, or orientation of each constituting microphase, i.e., neither sta- 
tistical homogeneity nor isotropy is assumed. Consider an arbitrary finite 
sample of volume Q and boundary dQ. The field variables for the problem 
are the stress, cr(x), and strain, e(x), tensors and the displacement vector, 
u(x). The constitutive equation relates the stress to the strain through the 
stiffness, C(x), or the compliance, D(x) = C^"*", tensors; in what follows, the 
X dependence is implicit, 

o- = C : £; £ = D : <T (1) 

Strain is related to displacement through the kinematic relation (field 
equation). 



4 



e = 1 [Vu + V^u] (2) 

For elastostatically admissible boundary data, effective properties (effec- 
tive stiffness and compliance tensors) are defined by relating the domain 
averages of the stress and the strain tensors, 

(o-)n = C^ff : {e)n = D^/^ : {ct)^ (3) 

where the domain average is defined by, 

(Q)n = ^ [ QdQ 

Effective properties defined above, depend upon the boundary conditions 
under which the strain and stress fields are generated. It is, therefore, of 
interest to study if there exist special boundary conditions which bound 
the effect ive properties associated with any b oundary data. It was shown 



in Refs. 



Nemat-Nasser and Horil (119951 . Il999h that the elastic energy and 



complementary energy associated with the elastostatic system are bound 
by special boundary conditions and that this fact could be used to place 
strict bounds on the effective parameters. To this end, the following two 
stress/strain states were defined: 

• Weakly Kinematically Admissible Strain Fields: Any compatible strain 
field with a prescribed average value. 

• Weakly Statically Admissible Stress Fields: Any compatible stress field 
with a prescribed average value. 
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And the following two universal theorems were proved: 

• Universal theorem for elastic strain energy: Among all weakly kine- 
matically admissible strain fields, the strain field produced by uniform 
boundary tractions renders the total strain energy an absolute mini- 
mum. 

• Universal theorem for complementary elastic energy: Among all weakly 
statically admissible stress fields, the stress field produced by linear 
displacement (uniform strain) boundary data renders the total comple- 
mentary energy an absolute minimum. 

In the present paper we show that analogous universal theorems exist for 
the elastodynamic case. These universal theorems may subsequently be used 
to place strict bounds on the effective dynamic parameters. 

3. Elastodynamic Universal Theorems 

Now consider the dynamic case of a general heterogeneous solid which 
consists of various elastic phases. As for the elastostatic case, there is no re- 
striction on the number, geometry, material, or orientation of each constitut- 
ing microphase, i.e., neither statistical homogeneity nor isotropy is assumed. 
Consider an arbitrary finite sample of volume Q and boundary dQ. The field 
variables are represented by the stress, cr(x), and strain, e(x), tensors and 
the vectors of momentum, p(x), and velocity, u(x). Constitutive relations 
relating the stress to the strain are given by Eqs. ([T]). The relation between 
the momentum and the velocity at every point in Q is. 
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p = pu (4) 

where p(x) is the density. The field equations relate the stress to the 
momentum and the strain to the velocity, 



V • o- = p 

£ = (Vu + 



(5) 



Dynamic homogenizat i on is an active area of research (Se e 



Norris and Santosal fll992[ ): lNorrisl fll993h : IWang and Rokhlinl fl2002h : 



Norris ( 



1992) 



Andrianov et al 



(120081 ) and references therein.) Effective dynamic parameters are defined by 
relating the domain averages of the field v ariables. The general form of the 



averaged constitutive relation is given bv (IMilton and Willis! (120071): 



f 2009 ): lNemat-Nasser and Srivastava (,20 111 ): IWillisI (120111 ) 

(boii), 



Wilhs 



Srivastava and Nemat-Nasser 



{s)n = D : + ■ (u)^ 



(6) 



All the effective constitutive parameters are non-local in space and time. 
The effective parameters may be complex even if there is no dissipation in the 
system, the imaginary parts resulting from the asymmetries of, e. g., the unit 
cell of a periodic composite. D is the fourth-order effective compliance tensor 
which has minor symmetries, Dij^i = Dji^i = Diji^. It does not have the 
major symmetry associated with the elasticity or the compliance tensor but 
has a hermitian relationship over the major transformation, Dijki = [Dkuj]*, 
where * denotes complex conjugation. Effective density is a second-order 
tensor with a hermitian relationship over the transformation of its indices. 



Pij — [pA*-! ^^'i third-order coupling tensors with the relationship 

Cl = [C2 1* 
ijfc I kiji ■ 

Defined as above, the effective parameters depend upon the boundary 
conditions under which the domain averages of the field variables have been 
calculated. It is, therefore, of interest to investigate the existence of special 
boundary conditions which may bound the effective parameters by bounding 
the associated total elastic energy plus the kinetic energy and the total com- 
plementary elastic energy plus the kinetic energy for any boundary data. We 
present two theorems which are the elastodynamic equivalent of the elasto- 
static universal theorems stated above. These theorems prove the existence 
of special boundary conditions under which the total elastic energy plus the 
kinetic energy and the total complementary energy plus the kinetic energy 
achieve their absolute minima. 

3.1. Universal Theorems 

We define the following concepts of weak admissibility, 

• Weakly Kinematically Admissible Strain Fields: Any compatible strain 
field with a prescribed average value. 

• Weakly Dynamically Admissible Momentum Fields: Any compatible 
momentum field with a prescribed average value. 

It is not required for these fields to satisfy any specific boundary data, 
only their averages are required to be equal to specified values. 
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3.2. Universal Theorem for Total Elastic Strain and Kinetic Energy 

Consider a weakly admissible strain field e(x) which satisfies (e)^^ = e*'. 
The total elastic strain and kinetic energy associated with the body is given 

by, 



n(e;e°) = ^(e:C:e)n + ^(v-pv)n (7) 



where the inner product is given by, 



(e : C : e)n = ^ / e : C : e*drt 
1 r 

(v • pv)q = — / V p\*dVt 



(8) 



For any other weakly kinematically admissible strain field, e(x); (e(x))n = 
e*^, the total energy is given by, 

n(6;6°) = i(e:C:e)n + ^(v-pv)n (9) 

where v(x) is the velocity field resulting from the displacement field u(x) 
which corresponds to the strain e(x). Since C has major symmetry, the 
difference in the total energy for the two weakly admissible strain fields is. 



n(e; e°) - n(e; 6°) = ((e - e) : C : (e - e))n + -((e -e):C:e)n 

1 1 
+ 4((v - v) • p(v - v))n + -((v - v) • pw)a + c.c. 

(10) 

where c.c. denotes the complex conjugate of the preceding expression. 
Setting & — C : e and e = (e — e), the second term on the right hand side 
can be written as. 
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^((e - e) : C : e)Q = ^ / e^a^dfi 



Now we have, 



For the static case the last term in the above equation is equal to zero 
based on the conservation law a*j j = 0. For the dynamic case, however, we 
have cr*j j = p*, hence. 



i((e-e):C:e)^ = -^ / V ■ (6-* ■ - / p* ■ udQ 

1 /■ ^ 1 /■ - ^ ^ 

= —— / n ■ CT* ■ uddfl — / p* ■ udfl 

2" JdQ 2iz 



Now the last term in Eq. ( fTOl) is. 



i((v-v)-pv)^ = i^ /■ {i.p*rf^] (13) 

This cancels with the volume integral in Eq. f|T2|) if the two cases for 
which the energies are being compared are in harmonic motion with the 
same frequency. Moreover if the boundary conditions associated with the 
strain field e are such that, on the boundary, 

n(x)-6-(x)|an = n(x)-S (14) 

where I] is a constant tensor, then the surface integral in Eq. f|T2l) can 
be written as, 

10 



i / n &*- uddQ = ((e - e))n : S* (15) 
" Jdn 

This is zero when both e and e are weakly kinematically admissible fields 
with volume averages equal to e°. Similar considerations apply to the com- 
plex conjugate parts of Eq. ( ITOl) and it can be shown that the remaining 
terms in Eq. (flOl) are always real and positive given the symmetric, positive- 
definiteness of C and the scalar nature of p. 

The above treatment shows that among all weakly kinematically admissi- 
ble strain fields, the sum of the elastic strain and the kinetic energy contained 
in a finite body, in harmonic motion at a common frequency, is minimum for 
the case when the boundary conditions are one of uniform traction (constant 
stress). To summarize, our first universal theorem for the total energy is. 

In elastodynamics, among all weakly kinematically admissible strain fields 
at a given frequency, the strain field produced by uniform boundary tractions 
renders the total strain energy plus the kinetic energy an absolute minimum. 

3.3. Universal Theorem for Total Complementary Energy 

Now consider a weakly admissible momentum field p(x) which satisfies 
(p)n = P°- The total complementary elastic energy plus the kinetic energy 
associated with the body is given by, 

n'=(p; P°) = ^(o- : D : (t)o + ^(v ■ pv)^ (16) 

where D(x) = [C(x)]~^ is the compliance. For any other weakly ad- 
missible momentum field, p(x); (p(x))f7 = p", the total energy is now given 

by, 



11 



mp-, P°) = ^(^ ^ D : + i(v ■ pir)n (17) 
The difference in the total complementary energy of these states is, 



(p; p°) - n^(p; p°) = h(a-&):B:ia- + - a) : B : &)n 

1 1 

- v) ■ p(v - v))n + -((v - v) ■ pv)f^ + c.c. 

(18) 

Denoting e = C : <t and = (cr — ct), the second term on the right hand 
side can be written as, 



\(a-&):B:a)n = ^ I a^.e^^dQ (19) 



n 



2'' ' ' 2n 

As shown in the previous subsection, we have 



{^iju*),i + aijju* (20) 



so that Eq. (fT9|) becomes. 



—— [ n • cr ■ u*ddQ — / p ■ u*dQ 



(21) 



The volume integral in the above equation cancels the last term on the 
right side of Eq. (ITSll if the two weakly admissible cases are harmonic mo- 
tion with the same frequency. Therefore, for harmonic motion of the same 
frequency, say, u, if the boundary condition for the second case is such that 

ii(x)|aj, = U = -zcjU (22) 
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then the surface integral in Eq. fl2ip can be written as, 



- / n ■ a- ■ u*dafi = — di, .(in = — Pidn (23) 

Invoking the harmonic nature of the motion, we can write the final term 
in the above equations as, 

u* r ~ u* f 

- / = 77 / V^dn (24) 

The above integral goes to zero when the momentum fields for the two 
cases have the same average value. Similar comments apply to the complex 
conjugate part of Eq. f lTSj) . In light of this, our second universal theorem is. 

In elastodynamics, among all weakly dynamically admissible momentum 
fields at a given frequency, the momentum field produced by uniform boundary 
velocities renders the total complementary energy plus the kinetic energy an 
absolute minimum. 



4. Conclusions 



The static energy theorems given in Ref. iNemat- Nasser and Horil (Il999f ) 
proved that the elastic energy and the complementary energy of a hetero- 
geneous solid corresponding to any boundary conditions, are bounded from 
below by that corresponding to special boundary conditions, i.e., constant 
stress and linear displacement, respectively. Thus, the energy theorems show 
that the effective properties are bounded by the effective properties defined 
for these special boundary data. This was then used to calculate strict and 
computable bounds on the static effective parameters. For the dynamic case 
we, have proved the existence of analogous theorems which bound the total 
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strain energy plus the kinetic energy, and the total complementary energy 
plus the kinetic energy. It is worth to note that the boundary condition which 
provides a bound for the strain energy in the static case, i.e., the constant 
boundary tractions, also provides a bound for the total strain energy plus 
the kinetic energy for the dynamic case. On the other hand, the boundary 
condition which provides a bound for the complementary energy in the static 
case, i. e., linear displacements (or constant strain) does not provide a bound 
for the total complementary energy plus the kinetic energy in the dynamic 
case. Instead now, for a common average momentum, it is the uniform ve- 
locity boundary data that provide the bound for the total complementary 
energy plus the kinetic energy of the elastic composite. The energy bounds 
proved for the elastodynamic case in this paper are expected to provide strict 
bounds on the effective dynamic properties of heterogeneous composites. 
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